a='/home/wangdy/MNVTFdb/Server/Services/Predict/tmp'
b='/home/wangdy/MNVTFdb/Server/Services/Predict'
/home/caow/miniconda3/envs/mnvtfdb/bin/python $b/scripts/deltasvm_mnv.py $a/$1/$1.ref.fa $a/$1/$1.alt.fa $b/resources/models.weights.txt $a/$1/$1.tsv $b/resources/threhsolds.pbs.tsv 
sort -k1,1 -k3,3 -o $a/$1/$1.tsv $a/$1/$1.tsv
grep -v None $a/$1/$1.tsv >$a/$1/$1.filter.tsv
/home/caow/miniconda3/envs/mnvtfdb/bin/python $b/scripts/filter_promoter.py $b/resources/hg19_refseq_NMgene_promoter.sorted.bed $a/$1/$1.filter.tsv $a/$1/$1.f_prom.tsv
sed -i 's/_/\t/g' $a/$1/$1.f_prom.tsv
cat $b/head  $a/$1/$1.f_prom.tsv  >$a/$1/$2
